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We study the behavior of ultracold bosonic gases in the critical region above the Bose-Einstein 
condensation in the presence of an artificial magnetic field, B a rt . We show that the condensate 
fluctuations above the critical temperature T c cause the fluctuational susceptibility, Xfl? °f a uni¬ 
form gas to have a stronger power-law divergence than in an analogous superconducting system. 

Measuring such a divergence opens new ways of exploring critical properties of ultracold gases and 
an opportunity of an accurate determination of T c . We describe a method of measuring Xfl which 
requires a constant gradient in B a rt and suggest a way of creating such a field in experiment. 
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Amongst intensive simulation of condensed-matter ef¬ 
fects in cold atomic gases (see m for reviews), con¬ 
siderable attention was focused both on similarities and 
on striking differences in properties of superconducting 
systems on the one hand and ultra-cold Bose systems on 
the other (see [4] for review). Yet, the impact of fluctua¬ 
tions of the condensate order parameter above a critical 
temperature T c remains to be observed in atomic gases. 

In the vicinity of T c , i.e. for \r\ <C 1 where r = T/T c — 1 
is a reduced temperature, superconductivity can be de¬ 
scribed within the Ginzburg-Landau mean-field (MF) 
theory [5]. Its tremendous success for conventional clean 
superconductors is based on irrelevancy of the fluctua¬ 
tions for all achievable temperatures due the smallness 
of the Ginzburg number, Gi rsj 10 12 4- 10 14 . Here the 
Ginzburg number Gi defines the temperature interval, 
|t| < Gi, where fluctuational effects dominate [6]. How¬ 
ever, Gi is much larger in dirty superconductors so that 
temperatures r ^ Gi become attainable. In the tem¬ 
perature interval Gi < r <C 1 the MF results still domi¬ 
nate but fluctuational corrections become observable and 
lead to a sharp power-law r-dependence of conductivity 
[7] and magnetic response [8] above T c . The observa¬ 
tions made in Refs. 0 i were in excellent agreement 
with perturbative predictions by Aslamazov and Larkin, 
Maki, and Thompson mm- 

No similar observations exist for gases of cold bosons 
where analogs of the magnetic susceptibility and conduc¬ 
tivity are not readily available for measurements. On the 
other hand, the Ginzburg number Gi > 1 for a typical 
dilute cold bosonic gas: although it is proportional to a 
small gas parameter, the numerical coefficient is large, 
see Eq. © below. This makes the order-parameter fluc¬ 
tuations above T c strong and their effects potentially ob¬ 
servable. 

In this Letter we analyze the fluctuational contribu¬ 
tion, Xfl, to the susceptibility of a cold bosonic cloud 
in an artificial magnetic field, T> art , and suggest how to 
measure it. Up to now experimental studies of proper¬ 
ties of the BEC phase transition were mostly aimed at 


the divergent correlation length mm- Studying exper¬ 
imentally the critical susceptibility would allow one to 
measure another critical exponent thus building a more 
comprehensive picture of the phase transition. 

We show that the dependence of Xfl on the reduced 
temperature r is much sharper than in superconductiv¬ 
ity for the gas in a uniform trap, like that in Ref. m- 
We argue that a realistic measuring scheme can be based 
on using field B art with a constant gradient in space, and 
suggest a setup for creating such a field, see Fig. [TJ An 
implementation of such a scheme would expand the re¬ 
search in rotating condensates El and artificial gauge 
fields in general USM], which was mainly focused on in¬ 
creasing the flux densities to reach exotic states of matter, 
such as the quantum Hall regime [15] . 

An artificial magnetic field B art is created by imprint¬ 
ing an angular momentum and thus rotation on a cloud 
of neutral atoms. Neglecting interatomic interactions in 
the dilute cloud above T c , the corresponding change in 
the free energy of the cloud of radius R containing N 
atoms is To = — \mR 2 NujQ, which is equivalent to the 
free energy of rotation with frequency oo^ = CJ# /2 where 
ou b = B art /m is the analogue of the cyclotron frequency 
m- This looks like a rigid body rotation since the ran¬ 
dom thermal motion of atoms is averaged out. The sus¬ 
ceptibility per particle in natural units, 


1 d 2 F 0 
N du 2 B 


= ^rnR 2 , 


( 1 ) 


is proportional to the average moment of inertia per par¬ 
ticle, which is reduced below T c as the Bose-condensed 
part of the cloud does not contribute to it. The con¬ 
densate fluctuations above T c result in such a reduction, 
which we parameterize as fluctuational corrections to the 
susceptibility: X = Xo + Xfl- 

A detection of such a change in x requires measuring 
the ratio of angular momentum to angular velocity with a 
high precision. We suggest a scheme that removes the ne¬ 
cessity of difficult (if at all possible) direct measurements 
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FIG. 1. The proposed setup of creating a nonuniform artificial 
field, .Bart • Top panel - the generalized A-scheme: torque is 
imprinted on the cloud by difference in angular momenta of 
three laser beams that couple sublevels of the ground state, 

|pi, 2 ) with exited state |e). Detunings of the beams coupling 
1 02 ) and |e) change across the cloud due to a space-dependent 
Zeeman shift. Bottom panel - varying detunings results in 
different parts of a bosonic cloud acquiring different torques, 
which provides a gradient in the torque equivalent to VB art . 

by creating an artificial magnetic field with a constant 
gradient in the direction of the field, Fig. [l] 

The suggested scheme is based on modifying the stan¬ 
dard A-scheme pmn]. There the electronic ground state 
of an atom is Zeeman-split into two close sublevels, \gi^) 
with energies E\ ~£ 2 coupled to a single excited state, |e), 
of energy £ e by two laser beams, with the second being 
detuned by frequency 5 from the resonance [e e — 62 )/h. 
Angular momenta £ 1,2 are imprinted (e.g., with holo¬ 
graphic masks as in Ref. PU) on each beam so that their 
cross-sections have a Gaussian-Laguerre form with the 
Rabi frequencies adiabatically dependent on the atomic 
position, fip 2 (t) = |fii, 2 |(p/ poY 1 ^ 2 e~ p 2 / 2 Po e ^i, 20 . We 
put |fii| « |fi 2 | = fi, and i 2 = £ > 0 while £\ = 0. Then the 
coefficients of the internal atomic wave function, \^) = 
\gi) + b 2 (t)e~ l£2t / h |g 2 ) + b e {t)e~ l£et ^ h |e), in 
the rotating wave approximation obey the equations 

ibi#(t) = fii, 2 &e(*) , ib e (t) = fij&i(t) + fi 2 ^2 (t) • (2) 

In a steady state regime, each atom finds itself in the 
dark state HO], M(r)) = (fii(r) |^i) — 0 2 (t*) \g 2 ))/(Sl\/ 2 ), 
which is not directly coupled to the laser fields. Two 
other internal states orthogonal to | d{r)) are separated 
by a large gap and become redundant. A one-component 
wavefunction describing motion in the laser fields obeys 


a one-particle Schrodinger equation m in a vector po¬ 
tential A art = iA(d|Vd), corresponding to the artificial 
magnetic field in ^-direction [21] 

B “ = |T>A "l = f |uwIff /W - (3> 

Function /(5/fi) describes the field sensitivity to detun¬ 
ing; it equals [1 + S 2 / 2fi 2 ] -1 for \5\ <C fi and 2fi 2 /£ 2 for 
\S\ fi while its exact form is not relevant. 

To create a gradient in B art , we suggest to modify the 
standard scheme by coupling | g 2 ) and \e) with two, in¬ 
stead of one, laser beams carrying different angular mo¬ 
menta, £ a ^ and detuned by S a ^ from the resonance, see 
top panel in Fig. The gradient arises from linearly 
varying the Zeeman split (with a weak real magnetic 
field gradient in the ^-direction) between |0i, 2 ) and thus 
the ratio 5 a /5b, resulting in a different angular momenta 
transfer to different cross-sections along the beams. 

A rigorous description of the modified scheme amounts 
to replacing fi 2 in Eq. § by fi 2 (t) = fi a e Z(5at + fi^e 2 ^, 
where fi a ^(r) have the Gaussian-Laguerre form charac¬ 
terized by £ a ^. One cannot eliminate time dependence 
from thus modified Eqs. ©• However, we can demon¬ 
strate without exactly solving Eqs. © that an almost 
constant gradient of B art emerges with a proper choice 
of the parameters, as illustrated in Fig. [2] Such a gradi¬ 
ent results from setting the detunings in such a way that 

(i) 5 a = 0 and \5b\ ~ fi in the middle of the cloud and 

(ii) similar transverse profiles for fi a and fi^ are selected 
by choosing i a > i^ > 1 (see Fig. [|. In the Supplemen¬ 
tal Material, we illustrate that the values chosen for this 
figure are optimal for making this constant gradient ro¬ 
bust with respect to small changes in the applied real 
magnetic field and laser frequencies. 

Now we describe a response of the bosonic cloud with 



FIG. 2. (color online) The artificial magnetic field B art along 
the light-propagation axis z. The field is measured in units 
of Bo = fr/Po] the detuning is chosen so that 5 a — 0 at z — 0, 
5b = 0 at z = 1 and 5 a — 5b = fi (z is in arbitrary units). Each 
line represents a different combination of angular momenta: 
solid (red) is for £ a = 2 and ib — 3; dash-dotted (green) is for 
£ a — 2 and ib — 4; dashed (blue) is for £ a — 3 and ib — 4. 
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FIG. 3. (color online) The radial dependence, B art (p) of 
Eq. §, for Qb = 0 and zero detuning at different values 
of i a \ solid (red): t a — 2; dashed (blue): £ a = 3; dash-dotted 
(green): i a = 4. 

interatomic interactions to T> art created in the above 
scheme. As the interactions are typically too weak to ex¬ 
cite atoms from the dark state | d) [22] . an effective field 
theory can be formulated in terms of the one-component 
complex field T corresponding to | d). The action in the 
vicinity of T c is obtained from the full imaginary time 
action by retaining only the zeroth Matsubara-frequency 
component, T = T(r,cj n = 0), as higher frequencies do 
not contribute to singularities in r [23] . Thus results in 
[231 EU the universal classical action 

S W = / J dr ( /^/r -M^I 2 + ||^I 4 ) • (4) 

Here Sg = g — 2 gn, g is the chemical potential, n is the 
particle density at the center of the trap, g = Airh 2 a/m , 
and a is the scattering length. We assume weak in¬ 
teractions, an 1 / 3 < 1. The coupling constant g is 
practically unaffected by the A-scheme, which results 
only in changing the action by the minimal substitution, 
V —)> V — iA art , to include T> art [2T] . 

Essentially, S[Sfr] in Eq. @ is the standard Ginzburg- 
Landau functional [25] . However, for weakly interacting 
bosons in a homogeneous trap a r-dependence of fluc¬ 
tuations is sharper than in the superconducting second 
order phase transition. Below T c the chemical potential 
is fully governed by the interaction: Sg = —(3 /p)gnr, 
where p = 2,1 for a 3D gas in a homogeneous/harmonic 
trap. On the contrary, Sg above t c m is essentially the 
same as for the ideal gas: 

r f —ciT c t , trapped gas; , , 

^ | —C 2 T c r 2 , uniform gas. ^ ' 

where c\ « 2.2 and C 2 ~ 1.2 m- The interaction correc¬ 
tion to this is parametrically small for the harmonically 
trapped gas where gnr/\Sg\ ~ an 1 / 3 <C 1, and numeri¬ 
cally small in the region of interest for the uniform gas 
where gnr/\Sg\ < gn/(T c Gi) ^ 0.1, see Eq. § below. 


To show this, we rescale r —>> Xr and T —>> 6T to 
make the coefficients attached to all the three terms in 
the GL functional equal 1. The fluctuational weight 
e ~ s W/ T c becomes e~^ s where S is dimensionless and 
k = (2^\Sg\/gT c )(h 2 /2m) 3 / 2 . The Ginzburg criterion 
for suppression of the fluctuations [6] becomes u > 1. 
Substituting Sg and T c , we find k = (r/Gi) p / 2 . Thus for 
the most interesting case of the 3D gas in a homogeneous 
trap (p = 2) the Ginzburg criterion can be written as 

1 \r\ > Gi « 30an 1 / 3 . (6) 

This coincides up to a numerical prefactor with the condi¬ 
tion earlier formulated in the cold-atom context mm- 

The gas parameter an 1 / 3 is not small enough in typ¬ 
ical dilute gases to overcome the prefactor in Eq. © 
For example, for a typical density of trapped Rb atoms, 
10 12 T 10 13 cra -3 Gt an 1 / 3 is a few hundredths and Gi is 
just under 1. However, fine tuning the scattering length 
near a Feshbach resonance would allow one to reduce Gi 
by at least an order of magnitude [29], thus making the 
window available for observations. 

We consider cold atoms trapped in an optical lattice 
forming a stack of Ni layers, where fluctuational effects 
are stronger than in the bulk. We assume that the laser 
beams and thus B art are normal to the layers. Then the 
minimal substitution affects only in-layer components of 
the gradient term in Eq. Q, 7lV|| -T fiVy — ^A art , while 
the normal component is replaced by J\^j — ff^+il 2 , 
where j enumerates layers separated by a distance d and 
J is a weak inter-layer coupling. After integrating e~ s 
over the fields T, one finds uni eh the fluctuation con¬ 
tribution to free energy at r<l as follows 

p _ $T C ^_ ttT c _ 

$0 ££ n S ^(T) + fkUB(n + \) + 4Jsin 2 (Y) ’ 

Here is the total flux of B art through the layer, <£>o = 
2tt h plays the role of the elementary flux, n = 0, 1 , 2, ... 
labels the effective Landau levels in the field H art and k z 
is a quasi-momentum in the normal direction. The corre¬ 
sponding susceptibility per particle, Xft = — jjd 2 F/8 uj 2 b , 
is found in the weak-field limit [30] similarly to that for 
superconductivity 

Xfl/Xo = -(2c p /3N ll )[r p (r p ^g p )]~^ . (7) 

For the gas in a homogeneous trap (p = 2) the r de¬ 
pendence much sharper than in superconductivity results 
from the r 2 -dependence of Sg. Here Nn = 7V||(po) = 
TTndpQ is the number of particles in a magnetized part of 
a single layer and rg^ /{ci^T c ) is the anisotropy pa¬ 

rameter. Since the inter-layer coupling J is independent 
of other parameters, both the regimes g < Gi or g > Gi 
are possible. In the latter case a crossover between 2D 
and 3D behavior {r > g or r < g) lies in the region of 
the MF applicability, Eq. © 
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The fluctuations susceptibility in Eq. 0 is negative. 
Since the classical susceptibility xo is proportional to the 
moment of inertia, a reduction of the overall suscepti¬ 
bility in the critical region above T c is a fluctuational 
precursor of the nonclassical rotational inertia below T c 
(the Hess-Fairbanks effect, [19]). Similarly to supercon¬ 
ductivity, where the reduced magnetic susceptibility is 
a fluctuational precursor of the Meissner effect, such a 
reduction reflects the divergence of the size of the fluctu¬ 
ational superfluid droplets at T c . A similar reduction of x 
due to onset of superfluidity has recently been proposed 
32] for measuring the superfluid fraction below T c . 

Although the prefactor in Eq. 0 is small, a very sharp 
r dependence, especially in the case of the uniform gas, 
XflOcr _ 2 atr ^>77 (which is much sharper than in super¬ 
conductivity), makes the fluctuational effects observable. 
Such a sharp r dependence should be even more pro¬ 
nounced outside region 0 . Le - for r < Gi, where analyt¬ 
ical expression 0 is no longer valid but one still expects 
a critical behavior of x* In this case the appropriate crit¬ 
ical exponent can be in principle calculated numerically, 
as in the case of the critical correlation length [33] where 
it turned out to be in an excellent agreement with the 
experiment m- 

The next step is expressing Xfl via observable quan¬ 
tities. To increase the weight of Xfl in Eq. 0, JVj| 
and thus the laser beams aperture po should be rela¬ 
tively small. On the other hand, to ensure the linear 
response regime w.r.t. B art it should be large enough, 
port 1 / 3 >> 1, i.e. Hujb ^ h 2 /mpo <C T c , or equivalently 
ujbPo (T c /m ) 1//2 ~ vt- Taking n « 2.3 x 10 13 cm -3 as 
in measurements of the critical correlation length of a Rb 
cloud [TT], one can choose po ~ 10 n -1 / 3 ~ 3/im, (corre¬ 
sponding to ujb ~ 30Hz) to satisfy both the conditions. 

In a steady state the central part of each layer rotates 
(after averaging out atomic thermal motion) with its 
own angular velocity u )o = ^b/ 2 , which linearly changes 
from layer to layer due to the gradient of B art . After 
switching B art off, it is necessary to allow some time for 
the equilibration within each layer, i.e. for redistributing 
the angular momentum from the central, “magnetized” 
part across the layer by thermal collisions between par¬ 
ticles. The collision time, r co i ~ (a 2 m>r) -1 , can be ex¬ 
pressed via T c and Gi, Eq. 0 , so that in uniform gas 
t coi « (h/T c )( 30/Gi) 2 . This is about 30s if the scattering 
length a is tuned so that Gi « 0 . 1 . So one should use the 
Feshbach resonance again to temporarily increase a in 
order to facilitate the angular momentum redistribution. 

Expressing the angular momentum of the central part 
of the layer in terms of ujb using Eqs. ([l]) and ([7]), we find 
the angular velocity of each layer proportional to the field 
B art in this layer as follows 


N\\(Po)x(Po) UB 
N\\(R)x(R) 2 


Po_ 

R A 


1 + 


Xfl(po/r)l 
Xo(po) \ 2 ’ 


where the fluctuational corrections are included only in 


x(po) as their relative contribution is much smaller in 
x(R). Noticing that Xfl is negligible at r ^ 1 and ujb is 
T-independent, we find that in each layer 

= 1 + Xfl ( Po , t)/xo(po) ( 8 ) 

To extract Xfl °ne needs to measure uj with a high preci¬ 
sion, while keeping ujb under control. Having a constant 
gradient in £> art , which makes each layer to rotate at dif¬ 
ferent frequencies, achieves precisely that. 

To detect and measure such a differential rotation one 
can apply a short laser pulse to make each layer elon¬ 
gated, as in the superfluid case [34]. This elongation is 
preserved if the dephasing time due to atomic collisions 
is much longer than the rotation period. Restoring the 
scattering length to the value corresponding to Gi « 0.1 
gives r co i ~ 30s that would preserve the shape for hun¬ 
dreds of rotations with ujb ~ 30Hz as above. Then after 
rotating for time to ~ 7r/AcJo, where Acoo is a difference 
in angular velocities of the two outer layers, the orien¬ 
tations of all layers will be uniformly distributed over 
all angles and the projection of the entire cloud along 
the symmetry axis will change from the elongated to the 
round one. The T-dependence to(r) can be found by re¬ 
peated measurements of optical density along the 2 axis 
at different temperatures in the vicinity of T c . Then, in¬ 
verting Eq. one expresses the fluctuations suscepti¬ 
bility in terms of directly measurable dephasing times as 
Xfl/Xo = 1 — to(r)/to(l), which should reveal the critical 
temperature dependence 0. An additional experimen¬ 
tal control may be achieved by measuring the revival time 
Nito when all the elongated layers are aligned again. 

In conclusion, we have shown that an impact of the 
order parameter fluctuations on properties of ultracold 
bosonic systems at the onset of the BEC could be exper¬ 
imentally accessible. In contrast to superconductivity, 
the fluctuational susceptibility is fully described [25] by 
the GL functional 0 if the Ginzburg number (J6| is small 
enough. Thus measuring Xfl will provide a new way of 
studying bosons critical behavior near the condensation 
transition, complementing recent studies of the critical 
correlation length num]. Although Xfl can be detected 
for an ultracold gas in a usual harmonic trap [30j, the 
most striking effect is expected for the gas in a uniform 
trap like the one recently implemented in Ref. [13], which 
is the most suitable platform for studying the BEC phase 
transition [12]. In this case we have found the unusually 
sharp critical dependence of Xfl on T— T c , Eq. 0 with the 
corresponding critical exponent 7 equal 2 for r > 7 . Such 
a sharp criticality near the transition can provide another 
way to accurately determine T c . Finally, we hope that 
creating artificial magnetic fields with a constant gradi¬ 
ent will find other applications in ultracold systems. 

We are grateful to Peter Kruger and Mike Gunn for 
useful discussions. One of us (I.V.L.) gratefully acknowl¬ 
edges support from the Leverhulme Trust via the Grant 
No. RPG-380. 
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1 The chemical potential above T c 

We want to outline the reasoning that leads to the results in Eq. (6). Above the transition, nearly all the particles 
are in the excited states. This number is, at least above the transition, practically temperature independent as 
the number of particles in the ground state is vanishingly small. For a density of states of p a (e) = C a e a ~ 1 one 
finds 


n th = j dep(e) - l) 1 = C a T d T (a) Li a (e*) . (1) 

We used the definition of the polylogarithm 

1 f°° v a ~ l 

Li “ (x) = f (a)L ^7 ~l dy ' (2) 

At the critical temperature T c , the chemical potential /x is 0. We have to demand that small changes in 
temperature AT = T c r change /i in such a way as to leave the total particle number invariant. We use that 
(T c + A T) a « T c (1 + ar) and the expansion of the polylogaritm[? ] 

Li Q (e x ) = T(1 - a) (-x) a_1 + ^ ^ k h k , 

k 

where £(r) = Y^=i rF the Riemann zeta function. This leads to 

_ 1 Li a (er) - Li a (l) 

ol Li a (1) 

Since [i is small, we can focus on the relevant orders in the polynomial and obtain 

-Li a (l)ar = U a (e*)=T(l-a)(^j + ((a - 1) (3) 

For the trapped gas a = 3 and close enough to criticality the term quadratic in ^ and logarithmic corrections 
can be neglected, one finds 


T c 


303 ) 

C(2) 


t = -C 1 T c t 


(4) 


Note: It appears as if the polylogarithm diverges for integer a in the Gamma term. However, this divergence 
is exactly cancelled by the divergence of the £(1) term. A careful limiting procedure shows that for integer k 

the limit lim Q _> fc+ i T(1 - a)(-x) a ~ l + C(a fc 7 - x k = ^ (Yln=i « “ ln (~ x )) • 

For o = 3/2 (free uniform gas) and close to criticality, the linear in ^ term can be neglected and one has 


/i 1 /3C(3/2)\ 2 , , 

T'~~» f— T 


1 



2 The susceptibility in the uniform system 

The following derivation is similar to the treatment of the subject in the context of superconductivity in [? ]. 
Starting from the action in Eq. (5) we modify it to contain the magnetic vector potential A = In addition, 

the anisotropic stacked system (with spacing d) is introduced by replacing the gradient in the z direction by 
(V Z T) 2 J |\k n +i — T n | 2 ,where n labels the layers. We work in the Gaussian approximation, so the quartic 
term proportional g is neglected. To diagonalize, the Landau level representation is introduced 

V(r) = '£*n,kAn(r\\)e ik ’ d , 

n.k z 


where ry is the position vector within one layer and (j) n is the wavefunction of the nth Landau level (which 
is degenerate with the factor AH/4> 0 , A being the sample area and 4> 0 — ^ is the flux quantum) and k z is 
evaluated in the first Brillouin zone. Each state has the energy T c C 2 T 2 + Hujb (n + 1/2) + 2 J cos ( k z d ). The 
parameter ujb — B/m is the cyclotron frequency, and C 2 is just the constants calculated in (5). As a reminder, 
the partition function Z = f Z}\kL>T*e _,s ^’^ 1 can be evaluated for a quadratic bosonic action by using the 
standard relation f OTOT*e - ^ A ^ = det A -1 . Taking the logarithm yields the formula for the free energy, 
F = —T c In Z, given by 


F{r,h) 


AB 

$0 


T c y> 


r Ca +4 h(n+ 1/2) + (1 - cos (k z l)) 

7r/c a 


( 6 ) 


where h is a reduced magnetic field 4 h — hB/mT c C 2 - Since the sum is divergent one should introduce an upper 
cutoff n c « T c /Hujb ~ 1 /h and for convenience it is useful to introduce the notation f = r 2 + | (1 — cos ( k z d )), 
where rj — is the anisotropy parameter. Exchanging the sum and logarithm and using [? ] 


E(*) 


lim 

n c —>-oo 


z(z + 1) (z + 2) • • • (z + n c 


1 )’ 


(7) 


one obtains the approximation 

+ In [r (1/2 + t/A h)} - In \n c \ n// 4 ^ 1 / 2 ] . (8) 

Expanding to second order in h and converting the sum over the momenta into an integral and keeping only 
the contribution ~ h 2 leaves 


E n c ln 


4 hc n 


ANT C r /l ldk z 1 h 2 
<o J-tt/i 2tt 3 f 


ANTh 2 f" dd ( 2 , V „ 

nKrl.&h + 2 (I 



Performing the integration 


J 0 ( r2 + 2 ^ _ C0S =( t2 ( t2 + V )) 


- 1/2 


With V = ANl this gives for the free energy density 

F(t) _ Th 2 


/ W = 


v 3tt/£ 2 ^/r 2 (r 2 + rj) 


= --B 2 , 


(9) 


( 10 ) 


( 11 ) 


From which the susceptibility in Eq. (8) of the paper can be read off. 

It is however helpful to rewrite everything in units of the energetic susceptibility per particle 

~ _ 1 d 2 F _ m 2 Xfi _ 1 VV 

x ~ NV d(hw B f~ n 2 ~ T c y r 2 (r 2 + v ) 

j- 2 

where we used that T c ~ y/fj describes the suppression of the critical temperature due to the anisotropy, l xy 
being the interparticle distance within the layer. The above formula for the susceptibility interpolates between 
the strongly anisotropic(^ single layer for r rf) and for the quasi three dimensional case (r g ). Note that 
the result differs from the similar result in superconductivity, namely here the divergence of the susceptibility 
is stronger. However, a system with an additional trapping in the x — y direction behaves different. In the 
center of the trap where fluctuations are not influenced by the finite trap size, the correct susceptibility can be 


2 



recovered by letting r 2 —>> r. However, this really holds only in the center of the trap and for not too small r. 
In the thermodynamic limit the fluctuational behaviour is suppressed, as we will demonstrate elsewhere. 

One could speculate whether the Ginzburg criterion in the anisotropic system might considerably deviate 
from Eq. (7) of the paper. However, a similar rescaling as for the isotropic uniform system results in 


Gi 


35 


mJd 2 


i/e 

an 1 / 3 . 


Because of the small exponent, the Ginzburg criterion deviates only slightly for all values of practical interest. 
At this point it should be mentioned that in a real system the confinement may influence the effective scattering 
length. 


3 Gauge potentials for multilevel systems 

The following outline follows the review by ? ]. 

In the conventional time-independent A scheme, as exemplified by eq. (1) in the main text, the Hamiltonian H 
can be written as 


H -(£ +y w) ® i+u - 

where V (r) is the trapping potential and U is the spatially dependent matrix, 

( 0 fit 0 \ 

u = -\ si 1 o n 2 , 

\ 0 n* 2 0 / 

that describes the coupling of the internal states |pi/ 2 ) , \e) via the light. The spatial dependence lies in the 
f^(r). The matrix U can be locally diagonalized by the states |d(r)), |6 ± (r)), with 

\d(r)) = ±(Sl 2 \g 1 )-Sl 1 \g 2 )), 

|i,±(t|> "7i(Tr l! ' ,) + Tr l! ' 2)±|e) )^ 

The energy scale is given by £1 = |fb| 2 + 10 2 1 2 and the descriptions dark state for | d) and bright states for 

| b±) are derived from the coupling properties of the states to the light and the excited state 

The crucial point is that, after the preparation into one of the eigenstates, the atom moves only slowly in 
the light field. Then the internal state follows adiabatically the evolution of the state, say | d) , and does not 
(appreciably) transit into any of the other two states, |d=). The full state can at any point be represented in the 
spatially dependent basis \d ), |=b) 


Considering that 


i'l'(r)) - Yi V> x (r) lx(r)> • 

X=d,± 

V(^ x |x)) = (V^ x )|x)+^xlVx), 


( 12 ) 

(13) 


one can determine the effective action of the momentum operator P. Projecting the Schrodinger equation onto 
the state |d(r)), one obtains the effective equation for the amplitude 'ipd 



(P - A art ) 2 
2 m 


+ V + W d 


^d • 


The particle evolution of ipd is described as the original equation of motion, but with an added vector potential 
A art = ih(d\Vd) and the additional scalar potential W = ti 2 |(6 + + b~\Vd)\ /2m. In the following we assume 
that the trapping is stronger than the scalar potential; or, in the case of a uniform trap, that an external 
potential cancels its effect. 
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4 Detuning of the artificial magnetic field 

Whereas the previous section outlined in very general terms how an artificial magentic field is created, we want 
to describe here the more generalized setup from Fig. (1) in the main text of paper from the point of view 
of the dressed states of the atom. There are three internal states of the atom, two quasi-degenerate ground- 
states \gi) , |< 72 ) and one excited state |e), each having respectively the energy Si,s 2 and £ e in the unperturbed 
Hamiltonian. The spatially dependent light field couples \gi) with |e) with Rabi frequency f7i(r)e zu;it , aq being 
the frequency of the laser, whereas the second laser couples |< 72 ) with |e). The second laser is supposed to be 
a superposition of two light fields I7 2 = f7 a (r)e zu;ot + Qb( r )e luJbt , with respective frequencies cu a and 0 J 5 . If we 
denote by cj e i = (£ e — £ 1 ) /h and cj e 2 = (£ e — £ 2 ) /fr the Bohr frequencies between the levels, we obtain from 
the Schrodinger equation the set of equations for the bi of the state \ij)) = b\ \gi) ± 62 I02) + b e |e) 

d ( h \ ( e,/h 0 \ / b 1 \ 

ih— b e = hi £ e /h %L e -ivat + ?h. e -iv b t 6 e ( 14 ) 

V h ) \ 0 q^ e i“at + ^ e iu b t e 2 //j J \b 2 J 

We are looking first at the case when there are only two lasers coupling the three levels, one can always find 
a time-independent frame of reference. By introducing a slow moving variable 61 , one can separate out the fast 
moving component of the bi = e~ l6it bi , the bi can be made to be time independent. The Si have to be chosen 
appropriately. If the lasers are in resonance, uq = a; e i, w 2 — u e2 , then only Si = €i/h allows for true stationarity. 
The evolution of the bi can be described as the action of an effective Hamiltonian In the resonant case it 
becomes 


/ 0 Sll 0 \ 

H eS = ■= fil 0 fi 2 (15) 

2 \ 0 n * 2 0 / 

which is equivalent to the formulation of eq. (1). 

With fi = ^/|fii| 2 + \n 2 | 2 one finds the three eigenstates 

jM) = t ( n 2 |(?i) + fil I g 2 )), j6±) = (fil l.'/i) ± fi |e) + fi 2 I// 2 )) • 

The dark state | d) has eigenenergy 0, whereas the two bright states | b±) have energy ±Ml/ 2. This is just a 
reproduction of the situation in section 3. 

Next we are trying to find the asymptotics for the artificial magnetic field for small and large values of S/Q,. 
If not all the laser are tuned to resonance, there is a certain freedom in choosing the individual Si to achieve 
this. The physics are however not affected by such a choice. 

For convenience we choose it to act as a perturbation to the Hamiltonian of the form 

/ -1 0 0 \ 

SV = (5 0 0 0 . 

\ 0 0 +1 / 


Here S is the small (with respect to Q) detuning. Using standard perturbation theory, the dark state is changed 
in leading order (for |fb| = |fi 2 | — ^/a/ 2) to 


= 1 + 


2 fi 2 J 


- 1/2 


S (fiifij + fi*fi 2 ) 
\/2fi 3 


(IM-I b-)) 


The resulting magnetic potential for a system prepared in the dark state that movies adiabatically in the Rabi 
fields is given as (see [? ] ) A art = i (d\Vd) . The light fields are supposed to be of the Gauss-Laguerre type 

f \ ^ (^) where p is the radial coordinate and </> the angular one. Especially, since we are interested 

in the magnetic field in the z direction we are looking for A ar t 5 </> = (d^d^d^ 

It becomes clear, from the form of the bright states, that |6 + ) = d<t>\ 6_), as the term in which the two 
differ, ±17 |e), is </> independent. Because both, | b±) perturb | d) with an opposite sign, their direct contribution 

to the vector potential of the ground state cancels and the only change to the vector field comes from the lower 

/ 2 \ — 1/2 2 

weight of the dark state, aka - p d p | d) s = (1 ± ^ J - p d p \d). Hence B aTt (S) « 2 q 2 +6 2 ^o,art, where H 0 , art is 

the magnetic field of the dark state on resonance. 
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In the opposite limit and for large detuning (S/il ^ 1) the unperturbed Hamiltonian is 


f -S 0 0 \ 

H 0 = h\ 0 S e 0 , 

\ 0 0 s ) 

where the eigenstates are just the original states \g \), \g 2 ) |e). Adding now a small coupling via the Rabi fields, 
one obtains perturbations of the form 



2fl 2 

(^e + S) 


- 1/2 


( |e) + 3 %; lsi> + 




Similar states appear for the perturbation of the other two states. The influence on the magnetic field can 
be easiest seen in the case where only fl 2 is angular dependent. Then the strength of the magnetic vector 
potential^ (e^d^e^ will be dominated by the prefactor of \g 2 ) state and the magnetic field will be of the order 

Tl 2 /S 2 or H art (r, S) ~ yp-H 0 , ar t(r), just as a naive extension of the weak detuning result would suggest. 

These findings, namely that B aTt (S) « 2 n*+6 2 for 8/Q. < 1; and H art (r, S) ~ ^H 0 , art (r) for 6/tl > 1, 
suggest that the magnetic field is an analytic function of S/Q and that the magnetic field can be written as 


5 ar t(r,J) = 5 0> art(r)/(Vfi) > 


as in eq, (4) of the paper. 

Now we are looking at the case when the third laser is turned on ^ 0, but is also far detuned from the other 
laser on the same transition 8 a — 5^ fh The Hamiltonian becomes time-dependent when fA ^ 0, ^ uj a . 

The time dependance cannot be eliminated by a rotation of the states, one will always have a system that 
oscillates with frequency A = S a — 5^. One can generalize the previous decomposition of the solution into 61 , & 2 
and b e . Apart from an overall phase evolution, each of the c$ should be a periodic function of time, with period 
given by T = ^ (Floquet theorem or Bloch theorem for a time dimension). Formally one can Fourier expand 
the Ci(t) to contain all the overtones, e.g. b\ \gi) —» b^e lAnt \gi,n). Quite generally, the state of the atom 

can be decomposed into a Floquet basis, \ip) = e l£t i c ^e lAnt |i), where n is an integer of the Floquet basis 
and i is the index of the standard set of \gi) , | g 2 ) and |e). For very large detuning, only the time-independent 
n — 0 subspace is relevant. However by slowly decreasing the detuning, the other Floquet sectors come into 
play. The stationary Schrodinger equation for the slow variables b™ is (given TL a in in resonance) 

(e + nA) b[ n) = 

(e + nA) b^ = tt a b^ + 

(s + nA) &<"> = nib[ n) + 0^4 n) + 0^4" +1) - 

The perturbation V = f \ |# 2 , n) (e, n — 1 | + H 3 |e, n) (g 2 ,n + l\ couples the different blocks indexed by 

n, each with (large) offset nA. Without the coupling, every block can be diagonalized into the same set of 
eigenstates just as the time-independent problem without detuning. The coupling between different blocks is 


V n , n +i = ^2 | i,n) (i,n\V\n+l,j) (n+ 1, j\ = 9,* b | e,n) (g 2 ,n+l\ , 

i,j 

Vn+i,n — ^6 |# 2 5 n + 1) (e,n| . 


It is notationally convenient to separate the two indices i.e. |i,n) | i) |n), as we are interested in the 

coupling between the n. If one starts with a single laser field Tl a in resonance and in a dark state (\d) \n = 0)) 
and introdces the second light field by adiabatically lowering the detuning A = oo to a finite (but large) A, one 
obtains to lowest order in perturbation theory 


k 1 ) 




_ ^2 
n 

0 

IT 


\ 

| n — 0) 

/ 



V 


\Shf 

n 

AO* 

~TT 

n;n a 


\ 

I n = -1). 


where we used that « Aj. The state obtains a small component that is oscillating but also carries the 

spatial information necessary for an artificial magnetic field. The magnetic vector potential to leading order in 
^can be obtained as before by calculating the overlap (cj) l \ V^ 1 ) (after normalization). The perturbation adds 
an additional term to the magnetic vector potential proportional to the angular momentum of the beam f \e llb< ^ 
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and it is the same contribution as in the previous case of a single far-detuned laser, since O 3 « 0. Therefore 
this amounts to effectively adding the magnetic fields created by the two Rabi fields individually. 

One can argue that this perturbative expansion breaks down for O/A « 1 . However, by choosing the angular 
momenta of the beams sufficiently different, one can make the value of O/A, for which the second derivative 
of the magnetic field vanishes and the gradient is approximately constant over an intervall of order O as small 
as wanted, thereby making sure that the perturbation theory stays valid at all times. If one of the fields is 
in resonance, then the detuning has two effects. On the one hand the magnetic field created by the first field 
weakens. On the other hand, the atom gets closer to resonance with the second field, which is stronger by the 

factor (j^ • In particular, the difference between 5 a and 5b can be chosen such that downward curvature from 
the detuning with the first laser is compensated by the opposite curvature of the second field, thereby creating 
a steady gradient. 

Simple plotting shows, that for Al = 1 , the plateau of the gradient forms at O/A = 0.7, for Al = 2 at 
O/A = 0.4 and for Al = 3 at A/O = 0.32. 



upper graph (blue) A/O « 2.0. In the middle graph (red) A/O « 2.5 and in the lower graph (green) A/O « 2.8. 
One can clearly see how the finetuning of A/O allows for the compensation of curvatures by creating a plateau 
of width ~ O. 
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